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Abstract 

Let r be geometric tree graph with m edges and consider the second 
order Sturm-Liouville operator C[u] = {—pu')' + qu acting on functions 
that are continuous on all of T, and twice continuously differentiable in 
the interior of each edge. The functions p and q are assumed uniformly 
continuous on each edge, and p strictly positive on T. The problem 
is to find a solution / : T — >• M to the problem C[f] — h with 2m 
additional conditions at the nodes of T. These node conditions include 
continuity at internal nodes, and jump conditions on the derivatives 
of / with respect to a positive measure p. Node conditions are given 
in the form of linear functionals l\ , . . . , li m acting on the space of 
admissible functions. A novel formula is given for the Green's function 
G : T x r — > K associated to this problem. Namely, the solution to 
the semi- homogenous problem £[/] = h, £i[f] = for i = 1, ... , 2m is 
given by f(x) = J r G(x, y)h(y) dp. 



1 Introduction 



The Sturm-Liouville differential operator 

c[f] = -(pf'y + qf 



on an interval, appears in the analysis of many different types of models in 
the natural sciences. The problem C[f] = h or C[f] = vqf together with 
appropriate boundary conditions arise when considering Kirchoff's law in 
electrical circuits, the balance of tension in a elastic string, or the steady 



state temperature in a heated rod (see for example Kreyszig 1999|, Hjortso 
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and Wolenski] [2009| , [Guenther and Lee] [1996]). A more complete review of 



the mathematical theory can be found in |Zettl| (2005 



The extension of operator ([I]) to the case of a domain composed of in- 
tervals arranged in a graph has received recent attention for the last thirty 
years (see for example Merkov| |1985| , Roth 1984 , Below [1988 . A complete 



bibliographical review with historical notes can be found in 



'okornyi and 



Borovskikh 12004 



1.1 Physical motivation 

The application example that follows serves as the motivation for the present 
study, and arises from a problem in stability of populations of organisms in 



river networks [Ramirez 2011.] 



When considering the dispersion of solutes or organisms in a river or 
stream, one might consider the following advection-diffusion model 



du 
dt 



D 



du 



dx 2 dx 



(2) 



Here u(x, t) is the concentration per unit length of the dispersed quantity at 
position x S [0,1], and time t > 0. The coefficients D and v are constant and 
strictly positive and denote the diffusivity and water velocity respectively 



(see Lutscher et al. 2006 for a justification of this model). 

Consider now a mathematical model for the dispersion of the same quan- 
tity in a collection of streams arranged in river network. One then might 
consider a domain in the form of a tree graph T where each stream in the 
network corresponds to an edge e of V, and the stream junctions and bound- 
ary points are the nodes of the graph. An edge e can be parametrized as the 
interval (0,/ e ), the point x = corresponding to the downstream node of 
edge e. We denote by T the union of all edges, and by T the set composed of 
r and the nodes of the graph. The most downstream point of the network 
(i.e. its outlet) is the root node of T and is denoted by 4>. Let u : T — > M be 
the longitudinal concentration, with u e denoting the restriction of u to edge 
e. Then u satisfies the following evolution equation 



du 
dt 



A[u], u(x,0) = f(x) 



(3) 



where the differential operator A is given on each edge by 



A[h] 



D P h" — v P h' . 



(4) 
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The operator A acts on functions that satisfy certain regularity conditions 
inside each edge, but more interestingly, one must also specify conditions 
at the nodes of the tree graph. For the application in river networks, an 
internal node n is located where edges ex, e 2 merge to form edge eo- Since 
all nodes are assumed to be oriented downstream, one can talk about the 
value of u e and u' e at node n by taking the appropriate one-sided limits. In 
the particular application of dispersion on T, one requires continuity of the 
concentration: 

u eo (n) = u ei (n) = u e2 (n), i, j = 1, . . . , k, (5) 

and a flux matching condition 

piu ei (n) + p 2 u e , 2 (n) = p Q u eo (n) (6) 

for some nonzero coefficients p e , e C T. The set of nodes that have a single 
incident edge is called the boundary of T, and dT is used for its notation. The 
phenomena of dispersion typically imposes Dirichlet (absorbing) boundary 
conditions at the root node, and Neumann (reflecting) condition at all other 
(upstream) nodes. Namely, 

u(</>) = 0, u'(n) = 0, nedr\{(f>}. (7) 

Consider the following integrating factors 

ap(x) 



?(*):= exp{- / ^dy}, q{x) 



D(x) 



(8) 



where the integral on the definition of p is taken along the unique path con- 
necting the root <p and the point x £ T. The functions v and D are defined 
on T by taking the values v e and D e on edge e respectively. Calculation of 
the resolvent of A, involves inverting the operator (a — A) for an arbitrary 
a > 0. For A given in Q, it follows that 

C:=^(a-A) (9) 

has the familiar form of a Sturm-Luiville operator on each edge, 

C[f}e = -(pfe)' + qfe. (10) 

One is then interested in solving the problem 

C[f] = h, fe£(T) (11) 
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where S(T) is the set of functions that are twice continuously differentiable 
inside each edge, satisfy internal node conditions ([5]) and ([6j); and boundary 
conditions ([7]). It can easily be shown that problem ( |11[ ) has a unique 
solution if and only if it is non-degenerate, that is the only solution to the 
homogenous problem C[f] = 0, / G Dom(£) is / = 0. 

Let Q be some right inverse mapping of C namely £[<?[/]] = / for all 



admissible /. Then Q is the Green's operator for problem problem (11). 
Moreover, it will be shown that: 



Theorem 1.1. If problem (11) is non- degenerate, it has a unique Green's 



function. Namely, there exists a function G : V x V — > M. such that the 



solution to (11), f = G[h], is f(x) = J r G(x, y)h(y) dy. 



The proof of the existence of the Green's function G defined in Theorem 
1.1 follows from standard arguments. Here, the proof is obtained by simply 
giving an explicit formula for G(x,y). Uniqueness of the Green's function 
follows from the hypothesis of non-degeneracy. 



In Pokornyi and Pryadiev [2004 the authors provide a proof of existence, 
uniqueness and a formula for G for general graphs. The goal here is to 
present a new formula, both simpler and less expensive to compute, for 
the Green's function in the case T is a tree graph. The techniques here are 
elementary and based on the classical Lagrange's method for Sturm-Liuoville 



problems (see for example Guenther and Lee |1996| ). 



The organization is as follows. The next section settles the notation 
and defines the class of Sturm-Liouville problems to be considered. Finally, 
section [3] is devoted to the construction and the formula for the the Green's 
function. 



2 Sturm-Liouville problems on tree graphs 
2.1 Tree graphs and functions 

By a tree graph we understand a finite collection of edges embedded in M 2 , 
joined with nodes and containing no loops. That is, for any two points 
x, y in the graph, there exists one single path through the graph joining 
them. We assume that each edge e of the graph allows a sufficiently smooth 
parametrization, contains no self-intersections, and is finite, therefore can be 
considered as the interval e = (0, l e ). The collection of all graphs is denoted 
by T. At each endpoint of an edge is located a node of T. The set of nodes is 
N(T) and boldface is used to denote individual nodes. The graph, including 
its nodes, is denoted as T := T U N(T). 
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Points in T are denoted by the pair (e, x) with < x < l e , or by single 
letters if specification of the edge is not crucial. If n is a node, let i(n) denote 
the set of incident edges at n, namely those for which n is an endpoint. 
Boundary nodes are those n with #i(n) = 1. The set of all boundary 
nodes of V is dT. The set of internal nodes is 1(F) = N(T) \ dT. Node n 
has #i(n) ^ 1 possible representations: for each e G i(n), one either has 
n = (e, 0) or n = (e,l e ). The representation of points in T is therefore 
dependent on the parametrization direction of its edges. 

The value of a function / : T — > R at a point in V is denoted as f e (x) = 
f(e, x). That is, f e is the restriction of / to the edge e. For a node n G dT 
located at the endpoint of some edge e, f(n) denotes the appropriate one- 
sided limit of f e . For an internal node n with i(n) = {e±, . . . , e n } the value 
f ei (n) denotes the one sided limit of f Ci as x approaches the endpoint of 
at which n is located, i = 1, . . . , n. If all these limits coincide, / is said to 
be continuous at n, and f(n) is defined as the common value. 

We must also differentiate functions given on T. For a point (e, x) £ T 
with < x < l e , the derivative f'(e,x) = f' e (x) is computed as the usual 
derivative of the restriction f e at x according to the particular parametriza- 
tion direction of e. A change in the orientation of the parametrization of the 
edge implies a sign change on f' e . Note that the sign of f" or {pf' e )' remains 
unchanged. For a node n located at an endpoint of edge e, we introduce 
the boundary derivative as the derivative "out of node n into edge e" : as 
if the parametrization of e has n = (e,0). Boundary derivatives are useful 
because they make the following equality hold 



regardless of whether the integral is computed from to l e , or from l e to 0. 

The space of functions that are n times continuously differentiable in T, 
is denoted by C n (T), n = 0,1,2,...; C{T) := C°{T). Clearly, such spaces 
are identifiable with direct sums of the form (J) e C n (0, Z e ). The set C(T) is 
composed of functions in C(T) that are also continuous at each node. 

2.2 Sturm-Liouville operators 

Let p,q €: C(T) be bounded with mfrp(x) > 0. The object of this study is 
the following differential operator 



where Vl(T) denotes the space of functions / G C(T) such that (pf)' G C(T). 




(12) 



C[f] :=-(pf')'+qf, f£V 2 p {T) 



(13) 
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Green's functions for Sturm-Liouville operators are useful for solving 
more general problems than the one outlined on the introduction. For this 



more general treatment we follow Pokornyi and Pryadiev [20041. 



If T contains m edges, then the dimension of N(C) := {u :G Vp(T); C[u] = 
0} is 2m. A basis {ipi, . . . , <f2m} for N(C) can be found as follows. Let e be 
the i-th edge and define </?2i-i and if2i as the solutions to —(p e f')' + Qef = 
on e satisfying 

PM-i(0) = 1, ^ 2i _i(0) = 0, ip 2i (0) = 0, ^(0) = 1 

extended to all of T via <^2i-i(e, x) = (f 2 i(e, x) = 0, < x < lg for all e^e. 
Consider now a collection of 2m linear functionals {£i, . . . ,^2m} defined on 
Pjj(r). The problem 

/€X»J(r), £[/] = /i, ^[/] = Ci, i = l...,2m (14) 
will be uniquely solvable if and only if the homogenous problem 

feV 2 p (T), C[f} = 0, £i[f} = 0,i = l...,2m (15) 
has no solution except the trivial solution / = 0. In this case we say that 



problem (14) is non-degenerate. 



Non-degeneracy can be characterized as follows. Let A be the matrix 
defined by Aij = £{[<Pj], i,j = l,...,2m. Non-degeneracy is therefore 



equivalent to det(A) ^ 0. In this case, the solution to problem (14) can be 



written explicitly. Let z be some solution to the semi-homogeneous problem 



zeV 2 p (T), £[z] = h, £i[z] =0, * = l...,2m 
then the solution f = z + Y^t=i a ifi to problem (14) must satisfy 

■ • • ¥2m 



(16) 



/ 

-0,2m 

Hence, Cramer's rule gives the useful formula 

<Pi 




z 
-ci 

~C2m 



(17) 



/ 



1 

det(A) 



det 



z 
-ci 

~C2m 



tP2m 



A 



(18) 
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2.3 The physical problem 

Motivated by physical applications, we now specialize to semi- homogenous 
Sturm-Liouville problems where the functionals : i = 1, . . . ,2m} corre- 
spond to a particular choice of conditions at the nodes of T. 
Consider the operator C[f] acting on the set 

£(T) = V 2 p (T) n C(T) n F p {T) n B D {T). (19) 

Where Bd(T) and F p (T) contain the boundary and weighted flux matching 
conditions respectively: 

B D = {/GC(f):/(n) = 0, near} (20) 
F P (T) = I f € C\T) : Pefe(n) = 0, n G I(T) \ . (21) 

I e6 (n) I 



The function p in (21) is assumed constant on edges and strictly positive. 
Other boundary conditions types than Dirichlet - like those in Q- can be 
considered without major changes to the arguments that follow. 

The conditions encoded in £(T) can be cast in terms of linear functionals: 
let n G I(T) with i(n) = {ei, . . . , e^}, and define the functionals 

L,i[f] = /e l+ »-/e», i = l,...,k~l, (22) 
k 

L,k[f] = E^/e!H- (23) 
i=l 

For a boundary node n located at the endpoint of edge e, define simply 

in[f] = fe(n). (24) 

Relabeling gives a collection of 2m functionals {£i, . . . ,4m}, such that the 
problem of finding / E £(T) satisfying £[/] = h can be written as the 



semi- homogenous problem (16) 



3 Construction of the Green's function 



The goal is to arrive at a formula for the solution to problem (16). The 
solution to the associated non-homogenous problem will then follow from 
(p). 
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Definition 3.1. A Green's function for operator £ is a function G : TxT — > 
M such that for all h G Ran(C), the solution to problem (16) is given by 

f(x) = J G(x,y)h(y)dy. (25) 

The operator Q : h \-t J T G(x,y)h(y) dy is called the Green's operator. 

The first step is elementary and consists on verifying properties of the 
Wronskian of functions on T. For f,g G C 1 (r) the Wronskian W[f, g] is 
defined on an edge e of T as 

W[f,g] e = f e g' e -g e f' e (26) 

Lemma 3.2. Let f,g, h be functions in C 1 (r). 

a. If f,g G TJp(T) then Lagrange's identity holds on each edge e, 

f e C[g] e -g e C[f] e = --^( Pe W[f,g] e ). (27) 

b. Iff,ge B D , then W[f,g] G B D (T). 

c. Iff,geC(f)nF p (T), then Y, eeiin)Pe W[f,gUn) = for all neI(T). 
Here, the derivatives in the definition Watn are replaced by boundary 
derivatives. 

d. hW[f,g]-fW[h,g}=gW{f,h}. 

e. If f,g G £*p(r) with C[f] = C[g] = on some edge, then pW[f, g] is 
constant there. 

Proof. Statement (a) follows from a simple calculation and (b) is obvious. 
For (c), it suffices to use continuity and change derivatives to boundary 
derivatives, 

^ Pe W[f e ,g e ](n) = f(n) ^ p e g' e b (n) + g(n) ^ p e f' e \n) = 0. 

(d) is obtained by rearranging terms. To prove (e), compute (p e W[f, g] e )' = 
fe(Pe9e)' ~ 9e(PefeY ', use {pefj = -Qefe and (p e g' e )' = -q e g e to finally get 
(PeW[f,g] e y = 0. □ 

For the following definition, and subsequent formulas, assume without 
loss of generality that the parametrization of T is such that for all n G dT, 
we have n = (e, l e ) where e is the edge n belongs to. 



S 



Definition 3.3. Refer to figure^ Let (e,x) G T, 

a. The two connected components ofT\{(e,x)} are denoted T(e,x) and 
A(e, x) respectively. The point (e, x) is adjoined as a boundary node 
to f(e, x) and A(e,x). By convention, F(e, x) is taken as the tree 
that contains the node (e,0). There is an edge denoted by e in both 
f (e, x) and A(e, x); it is parametrized as the intervals (0, x) and (x, l e ) 
respectively. 

b. For an edge e, f (e) := U x£(0M T(e,x), and A(e) := U a . e(0j/e) A(e, x). 

c. As in the case of the full tree, T(e,x) - without the bar - denotes the 
collection of points inside edges q/T(e, x). Similarly for A(e, x) ,T (e) , A(e) . 

d. £ (r(e)) is the set of functions f G V$(T(e)) DC(f (e)) n F p (T(e)) such 
that f(n) = for nG dF(e)\{(e,l e )}. 

e. Similarly, <?o(A(e)) is comprised of functions f G Vp(A(e)) PlC(A(e)) H 
F p (A(e)) such that f(n) = for n G 0A(e) \ {(e, 0)}. 



Lemma 3.4. Let e be a fixed edge. 



If problem (16) is non-degenerate, then 

on 



there exists solutions ^ r(e) G £ (T(e)), V' A{e) G £o(A(e)) to £V r(e) 
T(e), anrf Cip A ^ = on A(e). These functions can further be chosen so 
that pW[^ r ^,V A(e) ] = -1 



on e. 



Proof. If e contains no nodes in <9r, choose a node n' ^ 5r(e) and let e' be 
its edge. If e contains a node in dT, make n' equal to that node, and e' = e. 
Rearrange the basis {ipi, . . . ,(f2m} so that ipi and 922 are supported on e', 
and the functionals . . . ,^2m} so that i\[f] = f(n'). Let A be the matrix 



defined in section (2.2). By the nondegeneracy of problem (16), there exist 
a solution a = (a±, . . . , 02m) to Aa = where denotes the M. 2m vector 
that has one in the first coordinate, and zero elsewhere. The function tp = 
Y^i=i a ifi is a solution in T>p(T)nC(T)r\F p (T) to C[tp] = on all of T and such 
that that ip(n) = for n G d(T) \ {n'}. The restriction of ip to T(e) serves 
as the r equi red function ip r ( e \ A similar construction applies for ip A ( e \ By 
Lemma 



3.2 



is constant on e, and the desired normalization 
can be achieved if this constant is not zero. Assume on the contrary that 
on e. Since the Wronskian vanishes, there is k 7^ such 



that V£ (e) 



^e A(6) . 



The function / := V> r{e) lr( e ) + ^ A(e) lrw would then 



be a solution to the homogenous problem (15) violating the assumption of 
non-degeneracy. □ 
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Figure 1: Schematic representation of a tree T with six edges. For e and x as 
shown, the black sub-tree on the upper figure in T(e, x), the gray sub-tree is 
A(e, x). The middle and lower figures depict trees T(e) and A(e) respectively 
with a schematization of the functions and ip^ e \ 
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Remark 3.5. The computation of the tp^ e \ip r ^ can be performed quite 
inexpensively. For a boundary n ode n' = (e',l e /) G dT, the solution t/j 
constructed in the proof of lemma 3.4 can be restricted to define ip r ^ for 
all nodes e such that either e = e' or e' does not belong to T(e). Similarly 
it can be used to define ip^ e ^ for all nodes e such that e' does not belong to 
A(e). This implies that the linear system Aa = has to be solved only 
#dT times. 

The specific form of the Green's function can now be written. 



Theorem 3.6. Assume problem (16) is non- degenerate. The following func- 
tion is a Green's function for operator C, 



G(x,y) = — x 

Pe 



V r(e) (y)^ A(e) (^) 5 yer( e ,x) 



x G e. 



(28) 



Moreover, this function is unique in the class of continuous functions on T 
that are continuous with respect to the first variable. 

Proof. Let h G Ran(£), and / G S(T) a solution to C[f] = h. Fix an edge e, 
and x G e. Applying Lagrange's identity (27) for ip r ^ and / and integrating 
over r(e, x) with respect to the measure p gives 



T(e,x) 



- [PaW[^ ie) ,f]aPa 
aCF(e,x) 



where the sum on the right hand side is taken over all edges a of r(e, x). 
Parts (b) and (c) of lemma 3.2 ensure that all terms in the sum cancel except 
for the value at (e,x), 



^(e) hdp = -p e { x )p e W[^ e \ f] e (x). 



(29) 



T(e,x) 



Similarly, Lagrange's identity for ?/; A ( e ) and /, gives 

i, K ^hdp = p e {x)p e W[^ e \f] e (x). 



(30) 



A(e,x) 



Multiply equations (29) and (30) by i(j A ( e \x) and ip r ( e >(x) respectively, add 
the resulting equations, and apply part (d) of lemma 3.2 to the right hand 
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side of the result. Finally, since — 1 on e, 

^\x)^\y)h(y)dp{y) + 

^(x)^\y)h{y)dp(y) = f e (x)p e . 



T(e,x) 



(31) 



A(e,x) 



Since T is a disjoint union of T(e, x) and A(e, x), the function G(x, y) defined 
in (28) satisfies definition (3.1). Let h G C(T) arbitrary. It will be establised 
now that h G Ran(£) simply by showing that / := Qh solves C[f] = h. Write 



f{x) = ^ A{e \x) / ^hdp + ^ e \x) / ^ e) hdp 



r(e)\e 



A(e)xe 



+ ip A ^(x) / V r(e) /idp + ^ r ( e )(x) / ip A ^hdp. 



(32) 



Applying C to the first two terms in (32 ) gives zero since C[ip A ^] = £[ip A ^] 
0. A routine calculation finally shows that 



£[/] 



-W[V r(e) ,V A{e) ]e + ^[^ A(e) ] I" V (e Wp + £h/; r(e) ] t ifj A{e) h 

JO Jx 



dp 



which yields C[f] = h. Lastly, the non-degeneracy of problem (16) and 
the fact that Ran(£) = C(T), imply the uniqueness of G as stated in the 
theorem. □ 

Remark 3.7. The construction of the Green's function in Theorem\3.6\has 



one particular important advantage over the one proposed by Pokornyi and 



Pryadiev, 2004' '■ ^ n that work, G(x,y) is given as 

2m 

G(x, y) = H(x, y)-J2 ti[H(;y)Mx) 



(33) 



i=l 



where H(x, y) is equal to the Green's function of operator C on (0, l e ) if 
x, y £ e, and equal to zero whenever x and y belong to different edges. The 
functions rji are solutions to £[r]i\ = 0, = Note that this formula 

requires solving Aa = a total of 2m times to compute G(x, y) at single 
pair of points (e x ,x), (e y ,y) ofT. Via formula (28), one needs only the 
functions ijj A ( ex *) and x/j r ( e ^ and therefore, the system Aa = must be 
solved only twice. On the other hand, formula (33) has the advantage of 
using H, which is a diagonal fundamental solution to C[f] = h. 
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